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ABSTRACT 

The north east shell of SN 1006 is the most probable acceleration site of high 
energy electrons (up to ~ 100 TeV) with the Fermi acceleration mechanism at 
the shock front. We resolved non-thermal filaments from thermal emission in the 
shell with the excellent spatial resolution of Chandra. The thermal component 
is extended widely over about ~ 100 arcsec (about 1 pc at 1.8 kpc distance) 
in width, consistent with the shock width derived from the Sedov solution. The 
spectrum is fitted with a thin thermal plasma of kT = 0.24 keV in non-equilibrium 
ionization (NEI), typical for a young SNR. The non-thermal filaments are likely 
thin sheets with the scale widths of ~ 4 arcsec (0.04 pc) and ~ 20 arcsec (0.2 pc) 
at upstream and downstream, respectively. The spectra of the filaments are fitted 
with a power-law function of index 2.1-2.3, with no significant variation from 
position to position. In a standard diffusive shock acceleration (DSA) model, 
the extremely small scale length in upstream requires the magnetic field nearly 
perpendicular to the shock normal. The injection efficiency (77) from thermal to 
non-thermal electrons around the shock front is estimated to be ~ 1 x 10~^ under 
the assumption that the magnetic field in upstream is 10 fiG. In the filaments, 
the energy densities of the magnetic field and non-thermal electrons are similar 
to each other, and both are slightly smaller than that of thermal electrons, in 
the same order for each other. These results suggest that the acceleration occur 
in more compact region with larger efficiency than previous studies. 

Subject headings: acceleration of particles — supernova remnants: individual 
(SN 1006) — X-rays: ISM 



1. Introduction 



Since the discovery of cosmic rays (Hess 1912), the origin and acceleration mechanism 
up to 10^^'^ eV (the "knee" energy) have been long-standing problems. A breakthrough came 
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from the X-ray studies of SN 1006; Koyama et al. (1995) discovered synchrotron X-rays from 
the shells of this SNR, indicating the existence of extremely high energy electrons up to the 
knee energy produced by the first order Fermi acceleration. Further, Tanimori et al. (1998) 
confirmed the presence of high energy electrons with the detection of the TeV 7-rays, which 
are cosmic microwave photons up-scattered by high energy electrons (the inverse Compton 
process) in the north east shell (the NE shell) of SN 1006. The combined analysis of the 
synchrotron X-rays and inverse Compton TeV 7-rays nicely reproduces the observed flux and 
spectra, and predicts a rather weak magnetic field of 4-6 //G (Tanimori et al. 1998, 2001). 

Since these discoveries, detection of synchrotron X-rays and/or TeV 7-rays, from other 
shell-like SNRs has been accumulating: G347.3— 0.5 (Koyama et al. 1997; Slane et al. 1999; 
Muraishi et al. 2000; Enomoto et al. 2002), RCW 86 (Bamba et al. 2000; Borkowski et al. 
2001b), and G266.6— 1.2 (Slane et al. 2001). These discoveries provide good evidence for 
the cosmic ray acceleration at the shocked shell of SNRs. The mechanism of the cosmic ray 
acceleration has also been studied for a long time and the most plausible process is a diffusive 
shock acceleration (DSA) (Bell 1978; Blandford & Ostriker 1978; Drury 1983; Blandford & 
Eichler 1987; Jones & Ellison 1991; Malkov & Drury 2001). 

Apart from the globally successful picture of DSA, detailed but important processes, 
such as the injection, magnetic field configuration, and the refiection of accelerated particles, 
have not yet been well understood. The spatial distribution of accelerated particles respon- 
sible for the non-thermal X-rays, may provide key information on these unclear subjects. 
Previous observations, however, are limited in spatial resolution for a detailed study on the 
structure of shock acceleration process and injection efficiency. Although many observations 
and theoretical models are made for SN 1006, these problems are still open issue (Reynolds 
1998; Aharonian & Atoyan 1999; Vink et al. 2000; Ellison, Berezhko, & Baring 2000; Dyer 
et al. 2001; Allen, Petre, & Gotthelf 2001; Berezhko et al. 2002). 

In this paper, we report on the first results of the spectral and spatial studies on the 
thermal and non-thermal shock structure in the NE shell of SN 1006 with Chandra (§ 3). 
In § 4.1 and § 4.2, we discuss the spectral analyses and determine the scale widths of the 
structures for thermal and non-thermal electrons on the base of a simple DSA with shock 
parallel magnetic field. We also derive the injection efficiency (77) of non-thermal electrons 
from the thermal plasma near the shock front. Based on these results, we discuss possible 
implications on the DSA process in the NE shell of SN 1006. In this paper, we assume the 
distance of SN 1006 to be 1.8 kpc (Green 2001). 



-3- 



2. Observation 

We used the Chandra archival data of the ACIS on the NE shell of SN 1006 (Observation 
ID = 00732) observed on July 10-11, 2000 with the targeted position at (RA, DEC) = 
(15''03™51!6, — 41'^51™18!8). The satellite and instrument are described by Weisskopf et 
al. (1996) and Garmire et al. (2000), respectively. CCD chips 12, 13, SI, S2, S3, and S4 
were used with the pointing center on S3. Data acquisition from the ACIS was made in 
the Timed-Exposure Faint mode with the readout time of 3.24 s. The data reductions and 
analyses were made using the Chandra Interactive Analysis of Observations (CIAO) software 
version 2.2.1. Using the Level 2 processed events provided by the pipeline processing at the 
Chandra X-ray Center, we selected ASCA grades 0, 2, 3, 4, and 6, as the X-ray events. High 
energy electrons due to charged particles and hot and flickering pixels were removed. The 
effective exposure was ~ 68 ks for the observation. In this paper, we concentrated on the 
data of S3 (BI chip) because this chip has the best efficiency in soft X-rays required for the 
spectral analyses and its on-axis position provides the best point-spread function required 
for the spatial analysis. 

3. Analyses and Results 

3.1. Overall Image 

Figure 1 shows the true-color image from 12, 13, and S1-S4 for the NE shell of SN 1006. 
The image is contrasted in the 0.5-2.0 keV band (hereinafter, the softl band; red) and in the 
2.0-10.0 keV band (hereinafter, the hard band; blue) and binned to a resolution of 1 arcsec. 
The fine spatial resolution of Chandra unveils extremely narrow filaments in the hard band. 
They are running from north to south along the outer edge of the NE shell, parallel to the 
shock fronts observed by Ha emission line (Winkler & Long 1997). These filaments resemble 
the sheet-like structure of the shock simulated by Hester (1987). The softl band image, on 
the other hand, has a larger scale width similar to the ROSAT BRl image (Winkler & Long 
1997). Many clumpy sub-structures are also seen in this energy band. 

3.2. Inner Shell Region 

To resolve the thermal and non-thermal components, we made a spectrum from a bright 
clump found in the softl band image, which is located in the inner part of the NE shell ( "Inner 
region" with the dashed eUipse in Figure 2). The background region was selected from a 
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region out of the SNR, as is shown in Figure 2 with the dashed hues. 

The background-subtracted spectrum shown in Figure 3 has many emission hues. We 
hence determined the peak energies of the 5 brightest hues with a phenomenological model, 
a power-law continuum plus Gaussian lines. The most strong line structures are the peak at 
0.55 keV and the hump at 0.67 keV. These energies are nearly equal to the Ka and Lyman 
a lines of He- and H-likc oxygen, hence are attributable to highly ionized oxygen. Likewise, 
the other clear peaks at 0.87, L31, and 1.76 keV are most likely He-like Ka of Ne, Mg and Si, 
respectively. However, in detail all the observed line energies are systematically smaller than 
those of the relevant atomic data. These apparent energy shifts have been usually observed 
in a young SNR plasma in non-equilibrium ionization (NEI). The "energy shift" in this case 
is due to the different line ratio of many sub- levels and/or different ionization states. The 
oxygen Lyman a is isolated from the other hues of different ionization states, hence the NEI 
effect gives no energy shift. Still we see apparent down-shift of the observed line energy from 
that of the laboratory data. He-like Ka lines are complex of many fine structures with the 
split-energy of at most ~ 25 eV (for He-like silicon). Although the energy shift of He- like 
Ka lines due to NEI should be smaller than this split-energy, the observed energy shifts are 
systematically larger than the split-energy. We therefore regard that the apparent energy 
shifts are due mainly to energy calibration errors, hence fine-tuned the energy gain to reduce 
by 3.8%, the average shift of the 5 brightest lines. We then fitted the spectrum with a thin 
thermal plasma model in NEI calculated by Borkowski et al. (2001a). The abundances of C, 
N, O, Ne, Mg, Si, S, and Fe in the plasma were treated to be free parameters, whereas those of 
the other elements were fixed to the solar values (Anders & Grevesse 1989). The absorption 
column was calculated using the cross sections by Morrison & McCammon (1983) with the 
solar abundances. Since this NEI model exhibited systematic data excess at high energy 
above 2 keV, we added a power-law component and the fit improved dramatically. Figure 3 
and Table 1 show the best-fit models (dashed and solid lines for thermal and power-law 
components) and parameters, respectively. 

Instead of the phenomenological power-law model, we applied srcut in the XSPEC 
package as a more physical model. Details of the srcut model fitting are given in § A.l. The 
best-fit i/rouoff is 9.2 (8.6-10.3) xlO^^ Hz, with better reduced of 389.0/215 than that of 
the power-law model of 447.9/215 (see Table 1). 

As for the thermal components, we also tried the fitting with a plane shock model 
(XSPEC model vpshock) plus either a power-law or srcut and found no essential difference 
from the case of an NEI model. 

Although these simple models globally follow the data very well, all are rejected in the 
statistical point of view, leaving wavy residuals near the line structure as shown in Figure 3 
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(lower panel). This may be caused by improper response function in energy scale and/or 
in energy resolution. We assumed that the photons are uniformly distributed in flux and 
in temperature through the whole source region. This simple assumption may also partly 
responsible for the above systematic error, because in reality the source region is apparently 
clumpy (see Figure 1) and may have different temperatures, abundances, and/or ionization 
time scales. Since our principal aim of this paper is to examine the spatial structures and 
the spectra of non-thermal component, we do not examine in further detail for the thermal 
model. In the following analyses and discussion, we use the physical parameters cited in 
Table 1 as a good approximation. 

The spectrum of the inner region clump is softer than any other regions in the NE 
shell, which indicates that the contribution of the thermal component is the largest. Nev- 
ertheless the thermal photons are only 0.02% of the non-thermal ones if we limit the en- 
ergy band to 2.0-10.0 keV (the hard band); the non-thermal photons in the hard band are 
8.9x10"^ cuts s"-*^, while thermal photons are 2.0x10"^ cuts s~^. Therefore, in the following 
spatial analyses, we regard that all the photons in the hard band are non-thermal origin. 

As for the spatial analysis of the thermal emission, we use the hmited band of 0.4- 
0.8 keV (hereafter; the soft2 band) to optimize the signal-to-noise ratio, in which K-shell 
lines from Hc-likc oxygen (0.57 keV) contribute a most fraction of the X-ray emission (see 
Figure 3). Even in this thermal-optimized band, however, the count rates of the thermal 
and non-thermal emissions are comparable: thermal photons are 8.1x10"^ cuts s~^, while 
those of non-thermal are 5.4x10"^ cuts s~^. 



3.3. The Filaments 

The outer edge of the NE shell is outlined by several thin X-ray filaments. For the study 
of these filaments, we selected 6 rectangle regions in Figure 1, in which the filaments are 
straight and free from other structures like another filament and/or clumps. These regions 
(solid boxes) arc shown in Figure 2 with the designations of No. 1-6 from north to south. 
Since the SNR shell is moving (expanding) from the right to the left, we call the right and 
left side as downstream and upstream following the terminology of the shock phenomena. 

Figure 4 shows the intensity profile in the hard (2.0-10.0 keV: upper panel) and soft2 
(0.4-0.8 keV: lower panel) bands for each filament with the spatial resolution of 0.5 arcsec, 
where the horizontal axis (x-coordinate) runs from the east to west (upstream to down- 
stream) along the fine normal to the filaments. We see very fast decay in the downstream 
side and even faster rise in the upstream side. 
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To estimate the scale width, we define a simple empirical model as a function of position 
(x) for the profiles; 



where A and xo are the fiux and position at the filament peak, respectively. The scale widths 
are given by and for upstream and downstream, respectively (hereafter, "u" and "d" 
represent upstream and downstream, respectively). Since the scale width of the filaments is 
larger than the spatial resolution of Chandra (~ 0.5 arcsec = 1 bin in Figure 4), we ignore 
the effect of the point-spread function. 



As already noted in § 3.2, the hard band (2.0-10.0 keV) fiux is nearly pure non-thermal 
origin (see also the next paragraph). We therefore used the hard band profiles for the study 
of the non-thermal X-ray structures. The hard band profiles were fitted with a function 
/''(x) -|- C'' (hereafter, "h" represents the hard X-ray band), where is the background 
constant, which includes the cosmic and Galactic X-ray background and non X-ray events. 
The fittings were statistically accepted for all the filament profiles. The best-fit models and 
parameters are shown in Figure 4 with the solid lines and in Table 2, respectively. 

We then made the spectra of the filaments within the scale widths: (xq — w][ < x < 
Xq + ifjj) in Figure 4. The background spectra were made from the off-filamcnt downstream 
regions. All the background-subtracted spectra are featureless (no line structure) and extend 
to high energy side, which were fitted with an absorbed power-law model with the best-fit 
parameters given in Table 3. We thus confirm that the hard X-ray profiles represent those 
of non-thermal X-rays. 

To increase statistics, we summed all the data of 6 filaments (the combined-filament). 
The best-fit power-law model parameters for the spectrum of the combined-filament are 
listed in Table 3. The spectrum was also fitted with a srcut model (see § A.l). The best-fit 
i^roUoflf and the other parameters are also listed in Table 3. We note that the srcut model give 
slightly better x^/d.o.f. than that of a phenomenological power- law model (see Table 3). 

We further spatially divided the combined-filament and analyzed each spectrum. We 
however found no significant difference between the downstream and upstream, nor within 
the downstream; the photon index was nearly constant along the x-axis in the combined- 
filament. 




(1) 



3.3.1. Non-Thermal Structure 
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3.3.2. Thermal Structure 

To examine the structure of the thermal components, we used the soft2 band profiles 
(the lower panels of Figure 4). Contamination from the non-thermal photons would be very 
large even in this optimized band (see § 3.2). Therefore we calculated the fiux ratio of the 
non-thermal photons between in the hard and soft2 bands using the spectral parameters 
given in Table 3. Prom the flux ratio and the best-fit hard band profiles, we estimated the 
non-thermal contaminations as are given with the dashed lines in Figure 4. 

After subtracting these non-thermal contaminations, we fitted the soft2 band profiles 
with a model of f^{x) + C^, where is the background constant in the same sense as C^. 
Note that although we use "s" as the abbreviation of the soft2 band, it actually represents 
thermal X-rays. 

From Figure 4, we see that most of the photons at the filament peak is non-thermal 
origin, hence the statistics becomes poor to determine the position of the thermal peak (xq) 
independently. We thus fixed Xq to the best fit peak in the hard hand Xq. Also we set the 
upper bound of the fitting parameter w\ to be 900 arcsec, the same as the radius of SN 1006 
(Green 2001). The best-fit models and parameters for the thermal components {f^{x)) are 
shown with dotted lines in the lower panels of Figure 4 and in Table 2, respectively. 

3.3.3. Non-Thermal versus Thermal 

Figure 5(a) shows the relation of the scale widths between the downstream and upstream 
sides for each filament. Although there is a large scatter, is systematically smaller than 
in both the non-thermal and thermal emissions. 

Figure 5(b) shows the relation between and w^. We find that is significantly 
larger than w^, whereas w^^ and are comparable with each other. The mean values are 
= 4.3 arcsec = 0.04 pc, w'^ = 1.3 x 10^ arcsec = 1.1 pc, w;Jj = 5.1 arcsec = 0.04 pc, and 
= 18 arcsec = 0.2 pc. Note that the minimum value of w|j is only 0.98 arcsec = 0.01 pc 
(filament 4). Maybe since their wide scatter is due to that these are the projected values of 
the possible sheet-like with wavy and/or curved shape. 
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4. Discussion 

4.1. Thermal Plasma 

Although the best-fit NEI model in Table 1 is rejected statistically, it globally fits the 
thermal emission of SN 1006 as shown in Figure 3. The temperature {kT—0.24: keV) is 
similar to the results obtained by Vink et al. (2000), Dyer et al. (2001), and Allen, Petre, & 
Gotthelf (2001), but lower than that by Koyama et al. (1995). Since the spatial resolution 

of Chandra enables us to remove the non-thermal photons from the thermal emission more 
accurately than the ASCA case, the present results should give a more precise description 
on the thermal plasma. Like the previous observations (Koyama et al. 1995; Allen, Petre, & 
Gotthelf 2001), heavy elements, in particular iron, are overabundant, which implies that the 
X-ray emitting thermal plasma is dominated by the ejecta from type la SN. The fact that 
the thermal emission is enhanced at the inner shell region also suggests the ejecta origin. 

From the emission measure (E.M.) and assuming a uniform density plasma of a prorate 
shape with the 3-axis radii of 140, 120, and 120 arcsec, we estimate the density in the 
inner region to be 0.36 cm~^. Then from the best-fit ionization parameter, r is 2.6 x 10^° s 
= 8.3 X 10^ yr, roughly consistent with the age of SN 1006. 

Even in the soft2 band profiles, the thermal components are not prominent (see Fig- 
ure 4) , which prevents us from high quality study for the morphology of the thermal plasma. 
Nevertheless, we found that the profiles of thermal filaments are largely anti-symmetric. The 
scale width is very sharp and comparable to w^, whereas shows a relatively large scale 
width. Although couples to xq, which is frozen to the best-fit value in the hard band (see 
Figure 5(b)), we can say that the thermal shock has very sharp rise in upstream. 

The downstream scale width w^^ is comparable to the shock width derived from the 
Sedov solution of about 75 arcsec = 0.7 pc for SN 1006 having a 15 arcmin = 8 pc radius 
(Green 2001). Therefore, the thermal filaments may nicely trace the density profiles of the 
Sedov solution. 



4.2. Non-Thermal Filaments 

In this section, we interpret the scale width of hard band X-rays in upstream, w^. 

The hardest spectral regions, the non-thermal filaments, are the most probable site of 
the maximum energy acceleration. From the observation alone, we can not judge whether 
the filaments are strings or sheets in edge-on configuration. Hester (1987) suggested that 



-9- 



thin sheet-like shock fronts are seen as filaments on the edge of the SNR. The filaments 
seen in figures of Hester (1987) resemble the structure in SN 1006. Furthermore, filaments 
should be also seen in the inner part of the shell if they are strings, however it is not the 
observed case. Therefore, we assume that the filaments have sheet-like structure normal to 
the shock direction. The depth of the sheet is unclear but would be similar to or smaller 
than the length of the filament. We thus assume the depth of the sheet to be about 1 pc in 
the following discussion. 

The most likely scenario of cosmic ray acceleration at the SNR shock is diffusive shock 
acceleration (DSA). The predicted results from this model however are highly dependent on 
many parameters of the magnetic field, such as the angle between the magnetic field and 
shock normal direction, magnetic field strength and its fluctuation. Here, we investigate 
the observed profiles based on a DSA model, and estimate the physical quantities such as 
diffusion coefficient, the direction of magnetic field, maximum energy of accelerated electrons, 
and the injection efficiency. Although a realistic condition may be more complex as is 
discussed by Ellison, Berezhko, & Baring (2000) and Berezhko et al. (2002), we at first apply 
a simplest DSA model that the back reaction of accelerated particles is neglected. 



4-2.1. Diffusion coefficient in upstream 

The scale widths in the upstream side are largely scattered in the range of 0.98-11 arcsec 
(Table 2). Since this large scatter would be due to the projected effect of possible sheet-like 
structure with wavy and/or curved shape, real widths should be smaller than the observed 
(projected) width. To be conservative, we however adopted the mean value of 4.3 arcsec, or 
0.04 pc at 1.8 kpc distance in the following discussion. 

We use the results of srcut, the fit of the wide band spectra from X-ray to radio band. 
The high energy roll-off froUofr srcut is the consequent of either age-limit (acceleration dom- 
inant), synchrotron loss, or diffusive escape (see § A.l). The best-fit t'loUoff at the filaments 
is 2.6^Q'y X lO^'' Hz, which can be converted to l.llo'3 keV. These are consistent with the 
ASCA result of 3.0tol x 10^^ Hz or 1.241°;°^ keV by Dyer et al. (2001). 

Since most of the non-thermal X-ray photons are observed in the downstream region, 
the synchrotron radiation is mainly due to the downstream region. Using the best-fit z/roUoff? 
we constrain the maximum energy of electrons £^max and magnetic field in downstream 
from the eq.(A6); 

E^^^B',' = O.Srt'oZ [ergsG°-^]. (2) 
In the case of the strong shock, the magnetic field in downstream B^ can be related to that 
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in upstream as 

= (cos^ 0^ + sin^ 0^) ^ B^, (3) 

where r and are the compression ratio and the magnetic field angle to the shock normal 
direction in upstream. 

The diffusion coefficients in upstream (K^) is estimated from eq.(A9) {u^ = Us) as 
following; 

= 3.1 X 10^^ [cm^s-^], (4) 

where the shock speed Ug is assumed to be 2600 km (Laming et al. 1996). The diffusion 
coefficient in upstream is given from Using eq.(AlO), we can derive for the electrons of 

J. 1^ / 2n sin 6'u \ Emax 

i^u^gCul^cos ^u + Y^j^c, (5) 

where ^u(> 1) is the fiuctuation of the magnetic field. 

We can derive following equations from eq.(2), (3), (4), and (5), the maximum energy 
(-Emax) and magnetic field in upstream (B^) are, 

_i 

E^^ = 37eu^(cos^^u + ^) ' 

X (cos^ 0^ + sin^ 0^)~~^ [TeV] , (6) 



B^ = 4.0 X 10-^Cu cos^ + 



5cl I .^.2^ , sm^0^Y 



2 



X (cos^ 0^ + sin^ 0^y'^ [G] . (7) 

Since SN 1006 is a type la SN located at a high latitude of 460 pc above the Galactic 
plane (/ = 14.6 at 1.8 kpc distance), the interstellar magnetic field would be smaller than 
the typical value in the Galactic plane of the order of 10 fiG. We therefore conservatively 
assume that B^ = 10 /xG, which is slightly larger than the combined results of the X-ray and 
TeV Gamma-rays (Tanimori et al. 2001). Then eq.(7) becomes, 

el (^cos^ 0^ + ' (cos^ 0^ + sm' = 0.25. (8) 

For the parallel magnetic field (^u = 0°), becomes smaller than 1. This unreahs- 
tically small value of confiicts to a usual DSA which assume a parallel magnetic field 
configuration. 
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An alternative scenario is oblique magnetic field to the shock normal in this region 
{9-a > 0°). For simplicity, we neglect the back reaction of accelerated particle, then the 
compression ratio becomes 4 for strong shock limit. The allowed region of 9^ is given from 
the eq.(8) as; 

> 82°. (9) 
Thus the magnetic field in upstream near the filaments is almost perpendicular. 
Assuming — 90°, r = 4, and — 10 //G, -Bmax is given from eq.(6), 

^max = r-i^u ^ ■ 0.231H1 [TeV] 

= 37^7 [TeV]. (10) 

When particles are accelerated very efficiently, their back reactions to the shock can not 
be ignored hence r becomes larger than 4. Even if the compression ratio is the largest of 
r = 7 (Ellison, Berezhko, & Baring 2000; Berezhko et al. 2002), the allowed range is > 80°. 
Therefore, we can safely predict that the magnetic field in the NE shell of SN 1006 is nearly 
perpendicular to the shock normal. 

As for the downstream region, the observed spatial profile seems to be incompatible 
with the solution derived by Blandford & Ostriker (1978). The maximum electron energy 
-f^max would be determined by the balance of the time scales between the accelerating and 
the synchrotron coohng. These time scales may depend on the structure and the fluctuation 
of the magnetic field along the shock normal; both time scales become smaller with larger 
magnetic field, while the former becomes also smaller with larger fiuctuation of the magnetic 
field. 

The shock fiow may compresses and partly stretches the magnetic field in the radial 
direction, which may produce highly disordered magnetic field with small fraction of radial 
component, as discussed by Reynolds & Gilmore (1993) with the radio polarization data, 
who reported that only 15-20% of the magnetic-field energy in SN 1006 NE shell in radial 
polarization, and most of the magnetic field is disordered with the scale smaller than 0.2 pc. 

To determine the magnetic field in downstream, we thus need require more complicated 
processes such as the history of the shock propagation and the non-linear effect, and many 
assumptions, which is beyond this paper and would leave for a future study. 
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4-2.2. Injection Efficiency 

NT 

We define the injection efficiency rj = where n^^ and nj are the number densities 
of non-thermal and thermal electrons in the filaments. The depth of the filament (sheet-like) 
is assumed to be 1 pc (see § 4.2.1), with uniform electron density. 

The non-thermal electron fiux (energy) are estimated using the method given in Ap- 
pendix (§ A.l), where the non-thermal X-ray fiux and spectra of each filament are taken 
from Table 2 and Table 3. We adopt the minimum energy of non-thermal electrons (injec- 
tion energy) Ey^in to be 0.24 keV (the temperature of the thermal plasma; see Table 1). The 
magnetic field and the maximum energy are unknown parameters, however we can suggest 
that the magnetic field in downstream is significantly larger than that in upstream. In this 
section, we adopt the magnetic field B to be perpendicular to the shock normal, — AB^ — 
40 /xG arbitrary and the maximum energy i?max to be 37 TeV from the discussion in § 4.2.1. 
The derived number density (n^^), the total number {N^"^), and the total energy {E^'^) of 
non-thermal electrons in each filament are summarized in Table 4. 

For the estimation of thermal electron flux (energy), we adopted the projected profile of 
the thermal X-rays (the soft2 band) given in Table 2 and the spectral parameters of thermal 
plasma given in Table 3. The resultant number density {n]:), total number (Nj), and total 
energy (Ej) of thermal electrons in each filament are given in Table 4. The injection efficiency 
7] (= is then obtained as is shown in Table 4. All the values rj are nearly identical in 
each filament of ~ 1 x 10~^ but are about 2 times larger than 5 x 10~^ derived from the 
ASCA data by Allen, Petre, & Gotthelf (2001), although they also assumed as = 40 /iG 
arbitrary. For comparison, we estimate rj with larger magnetic field of B^ — 40 /iG, and find 
a larger 77 than Allen, Petre, & Gotthelf (2001). They estimated the number density of non- 
thermal electrons from larger regions than the filaments, because ASCA could not resolve the 
filaments. On the other hand, we found that the non-thermal electrons are confined in the 
thin filaments, which are the sites of ongoing acceleration of the non-thermal electrons. We 
thus regard that the present Chandra result suggests that the injection occurs more locally 
and more efficiently and as a result, it must be more realistic estimation of rj than that by 
Allen, Petre, & Gotthelf (2001). 

The energy densities of the magnetic field, the thermal plasma, and the non-thermal 
electrons in the filaments are 6.4 x 10^^^ ergs cm"^, 2.6 x 10^^° ergs cm^^, and 6.9 x 
10^^^ ergs cm^'^, respectively. Thus, at the shocked region, the magnetic field and non- 
thermal electrons are in energy equi-partition and are slightly smaller than the thermal 
energy (about 30%). 

As suggested by Ellison, Berezhko, & Baring (2000), the non-thermal protons should 
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carry larger energy than electrons, hence particle energy becomes larger than that of magnetic 
and possibly that of thermal. Wc therefore must consider non-linear effects suggested by 
Ellison, Berezhko, & Baring (2000) and by Berezhko et al. (2002), however more quantitative 
scenario including the non-linear effects must be a future work. 

5. Summary 

(1) X-ray emissions from the NE shell of SN 1006 are found to be composites of filaments, 
clumps and more extended diffuse emissions. We resolved the thermal and non-thermal X- 
rays using the different spectral shape and morphology. 

(2) The spectrum of the thermal component has kT — 0.24 keV temperature and rela- 
tively small ionization time scale rieT = 1.1 x lO^'^ s cm~^. The chemical compositions are 
overabundant, especially in iron, which suggests that the X-ray emitting plasma originates 
from ejecta and the progenitor is type la. 

(3) The non-thermal components can be described with a power-law function with the 
photon index 2.1-2.3 in the narrow filaments and 2.5 at the inner region of the shell. 

(4) The structure of the filaments shows different characteristics in thermal and non- 
thermal X-rays. The thermal plasma has the scale width of 1 pc in downstream, similar to the 
shock width derived from Sedov equations. The non-thermal filaments show extremely small 
scale width of ~ 0.04 pc and w\ ~ 0.2 pc, in upstream and downstream, respectively. 

(5) In a diffusive shock acceleration model, the observed thin filaments requires a nearly 
perpendicular magnetic field with the angle between the magnetic field in upstream (assumed 
as 10 jiG) and the shock normal of larger than ~ 80°. The maximum energy of the electrons 
are 30-40 TeV. 

(6) The injection efficiency is estimated to be 77 = 1x10"^, suggesting that thermal 
particles are injected locally and very effectively. Then the energy density of non-thermal 
electrons becomes comparable to that of the magnetic field and about 30% of the thermal 
energy density. Thus non-linear effect of the shock structure and acceleration mechanism 
must be considered. 

We thank the anonymous referee for his/her helpful comments. Our particular thanks 
are due to M. Hoshino, T. Terasawa, T. Yoshida, and S. Inutsuka for their fruitful discussions 
and comments. A.B. and M.U. are supported by JSPS Research Fellowship for Young 
Scientists. 
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A. Appendix 

In this appendix, we briefly introduce the relevant software tools and equations, which 
are used for the discussions in the text. 



A.l. Electron Spectra and srcut Model 
The spectrum of non-thermal electrons accelerated by the diffusive shock is (Bell 1978); 



J NT 

^ = /.(£;^^ + mec2)(^^^V2mec2^^^)-^xexp(-£;^V£;_), (Al) 



where n^^ and E^^ are the number density and the energy of non-thermal electrons, re- 
spectively. 

The synchrotron radiation power per unit frequency from a single electron of energy E 
in a magnetic field B is 

P i/, a) = ^ F - , A2 

where a and F{x) are the pitch angle and the function given by Rybicki & Lightman (1979). 
The peak frequency Uc is 

'icE'^QB svaa , . 

^ 47r(mec2)3 ■ 

Convolving the eq.(Al) with (A2), we obtain the spectrum of synchrotron radiation in 
the pitch angle a as, 

roo J NT 

/(^,«)=/ d£;^^^P(i.,a). (A4) 

Averaged over the pitch angle, we finally obtain the synchrotron radiation energy per unit 
volume, frequency, and time, 

1 r 

/(i^) = 2 y dasina/(z/,a). (A5) 

The observed spectrum is fitted with this model spectrum using the Chandra software 
srcut wiih 3 free parameters (Reynolds 1998; Reynolds & Keohane 1999). The normalization 
constant and spectral index are determined so as to reproduce the flux and slope in the radio 
band. These parameters are converted to k and p in eq.(Al). The other fitting parameter 
i^roUoflf is a function of E'max in eq.(Al) and magnetic field {B) as is shown below (Reynolds 
& Keohane 1999); 
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i^roiioff = 0.5 X 10^6 ( ( _3l^] [Hz]. (A6) 
This equation gives constraint on the maximum electron energy E'max and magnetic field B. 

Since available radio data of SN 1006 have spatial resolution far larger than the scale of 
the X-ray filaments, we have no accurate radio flux nor index from the position of the X-ray 
fllaments. We therefore fixed the radio index to the poor resolution radio result oi p = 2.14 
by Allen, Petre, & Gotthelf (2001) and the radio fiux is treated as a free parameter for the 
present srcut fitting. Allowing the radio index to vary from p = 2.0 top = 2.2, we see no 
significant difference in the best-fit parameters within the statistical errors. 



A. 2. Diffusion versus Advection 

The spatial structure of the rclativistic electrons produced by the diffusive shock accel- 
eration across the shock is determined by the competing process of diffusion and advection. 
The advection time scale (rad) is given by 

Tad = -, (A7) 

u 

where w is the scale width of the spatial distribution of the relativistic electrons and u is the 
flow speed. 

The diffusion time scale (rdif) is given from the random-walk theory as; 

Tdif = ^, (A8) 

where K is the diffusion coefficient. 

In order that the high energy electrons are accelerated at shock front, electrons in 
upstream should be diffused back to the shock front against the advection to the downstream 

2 

side (shock flow). Therefore, Tad should be nearly equal to Tdif, hence ^ ~ We thus 
obtain; 

K 

Wu ~ — . (A9) 

For more exact formalisms, see e.g. Bell (1978), Blandford & Ostriker (1978), and Drury 
(1983). 

Let the mean free path parallel to the magnetic field be a constant factor ^ times the 
gyro radius Tg. Then, the effective diffusion coefficient along the shock normal is given as 
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(Jokipii 1987; Skilling 1975) 

K = ler,c(cos^^„ + f^), (AlO) 

= % (All) 

where e, Ef., B, and are the electric charge, the energy of relativistic electrons, and the 
magnetic field, and the angle between the magnetic field in upstream and the shock normal, 
respectively. 
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Fig. 1. — The true-color images of SN 1006 NE shell binned with 1 arcsec scale. Red and 
blue are 0.5-2.0 keV and 2.0-10.0 keV, respectively, both in logarithmic scale. 
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Fig. 2. — The close-up view of the 0.5-10.0 keV band image of S3 chip with J2000 coordinates, 
binned with 1 arcsec scale. The gray scale (the left bar) is given logarithmically ranging from 
1 X 10^® to 1 X 10^^ cnts s~^ arcsec"^ . The inner and background regions for the spectral 
analyses and the filament regions for the spatial analyses (No. 1-6) are shown with dashed 
and solid lines, respectively. 
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Fig. 3. — The upper panel: The background-subtracted spectrum of the inner region 
(crosses). Dashed hne and sohd hens are the best-fit thin thermal and power-law mod- 
els, respectively. The lower panel: The data residuals from the best-fit two-components 
model. 




10 



J4) 




y 

, 1 Mil,-/ 







160 10 



1 

'imJ 




h. i,iiiiili,iilk,i, „ 






r* 'i ii'Hi'i 1 ' 


— ^ 

J 


1 — 


— ' — 1 — ' — 1— 







-(3) 


i^i„,„„i ; 















30 40 





„,„ : 


lor ' " "•ffff 




v+^^i ; 







10 20 30 40 



40 eo 



10 20 30 40 



Fig. 4. — The profiles of the filaments in SN 1006 NE shell. Upper panels show the profiles in 
the hard (2.0-10.0 keV) band, whereas the lower panels in the soft2 (0.4-0.8 keV) band with 
the best-fit models (solid lines). The dashed lines in the lower panels represent non-thermal 
photons extrapolated from the hard band flux of the power- law (see the upper panels) . The 
dotted lines are the thermal component after subtracting the non-thermal contamination 
(dashed hues). Upstream is to the left and downstream is to the right. 
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Fig. 5. — (a) The relation between and w^. The close and open triangles are for the 
thermal and non-thermal emissions, respectively, (b) The relation between and w^. The 
close and open boxes are in the upstream and downstream, respectively. 



Table 1. The best-fit parameters of the spectrum of the inner regionf 



Parameters Best-fit value 



Power-law Model 

Photon Index 2.51 (2.48-2.53) 

Flux'^crgs cm"2 g-i] _ 3 3 ^ IQ-^^ 
NEl model 



Temperature [keV] 


0.24 (0.21-0.26) 


Abundances'^ 




C 


(<0.1) 


N 


(<0.03) 





3.3 (3.0-3.5) 


Ne 


4.8 (4.4-5.2) 


Mg 


51 (42-61) 


Si 


131 (121-140) 


S 


10 (6.8-13) 


Fe 


37 (28-46) 


riT [xlO® s cm~3] 


10.8 (9.9-11.1) 


E.A^.'^xlO-'^-'^cm-''*] 


2.1 (1.9-2.2) 


Flux'^ergs cm^'^ s^-"-] . . 


3.3x10-12 


iVn [lO^o H cm-2] 


9.0 (8.5-9.4) 


Reduced [x^/d-o.f] 


447.9/215 



''Parentheses indicate single-parameter 90% con- 
fidence regions. 

•"In the 0.3-10.0 keV band. 

'^Abundance ratio relative to the solar value (An- 
ders & Grevesse 1989). 

'^E.M. = rigV, where Ug and V are the electron 
density and the volume, respectively. 
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Table 2. The best-fit parameters of the profiles of the filaments.^ 



No. 


Ah 






Reduced 


A" 






Reduced 




[cnts arcsec"^]'' 


[arcsec] 


[arcsec] 


[xVd-o.f.] 


[cnts circsec"^]'' 


[arcsec] 


[arcsec] 


[xVd.o.f.] 


1 ... 


53 


1.9 


36 


119.1/105 


55 


1.7 


1.6 x 10^ 


104.0/94 




(49-58) 


(1.3-3.0) 


(30-44) 




(49-61) 


(1.3-2.2) 


(>92) 




2... 


63 


11 


22 


347.3/296 


1.0x10^ 


14 


93 


481.8/298 




(58-68) 


(8.7-14) 


(19-25) 




(95-1.2x10^) 


(9.8-21) 


(73-1.3 x 10^) 




3... 


66 


10 


5.5 


141.3/114 


44 


2.3 


1.5 x 10^ 


119.8/95 




(58-74) 


(8.2-13) 


(4.1-7.2) 




(37-52) 


(0.78-4.3) 


(>56) 




4... 


68 


0.98 


18 


103.2/99 


38 


0.12 


1.5 X 10^ 


106.1/81 




(62-74) 


(0.73-1.3) 


(16-21) 




(34-42) 


(not determined) 


(>78) 




5... 


75 


3.8 


8.0 


218.3/166 


33 


4.5 


62 


236.5/167 




(66-85) 


(2.8-4.9) 


(6.5-9.9) 




(24-41) 


(2.2-8.7) 


(43-1.0 X 10^) 




6... 


87 


3.5 


19 


99.9/86 


l.lxlO^ 


3.4 


1.3 X 10^ 


118.4/87 




(79-96) 


(2.1-5.4) 


(15-26) 




(1.0x10^-1.3x10^) 


(2.4-5.0) 


(86-2.5 X 10^) 





^Parentheses indicate single-parameter 90% confidence regions. 

the 2.0-10.0 keV band. 
t^In the 0.4-0.8 keV band. 



Table 3. The best-fit parameters of the spectral fittings for the filaments.^ 



Parameters 


1 


2 


3 


4 


5 


6 


Total 


Power-low model 
















r 


2.0 (1.8-2.2) 


2.4 (2.3-2.6) 


2.3 (2.1-2.4) 


2.3 (2.1-2.4) 


2.3 (2.1-2.5) 


2.3 (2.1-2.4) 


2.31 (2.29-2.33) 


ATh [1021 H cm-2] .... 


1.5 (0.9-2.1) 


1.8 (1.5-2.0) 


1.2 (1.0-1.5) 


1.6 (1.3-2.0) 


1.4 (1.1-1.8) 


1.4 (1.1-1.7) 


1.6 (1.6-1.7) 


Flux'^ergs cm~'^ s~^] . . 


3.7 X 10-13 


5.4 X 10-13 


3.7 X 10-13 


2.5 X 10-13 


2.3 X 10-13 


4.2 X 10-13 


1.8x10-12 


Reduced [xVd.o.f] . 


44.1/43 


232.5/164 


196.3/140 


139.8/125 


143.8/125 


173.4/130 


273.3/211 


srcut model 
















frolloff [xlO^'^Hz] 


13 (3.6-94) 


1.7 (1.1-2.3) 


2.8 (1.8-4.7) 


2.9 (1.5-6.2) 


2.4 (1.3-5.0) 


2.8 (1.7-5.0) 


2.6 (1.9-3.3) 


A^H [10^1 H cm-2] .... 


1.3 (1.0-1.7) 


1.4 (1.2-1.5) 


1.0 (0.8-1.1) 


1.3 (1.1-1.5) 


1.1 (0.9-1.4) 


1.1 (0.9-1.3) 


1.3 (1.2-1.5) 


Flux'^ergs cm~-^ s~^] . . 


3.6 X 10"" 


5.3 X 10-" 


3.6 X 10-13 


2.4 X 10-13 


2.2 X 10-13 


4.1 X 10-13 


1.8 X 10-12 


Reduced [x^/d.o.f] . 


43.7/43 


225.1/164 


190.9/140 


139.1/125 


141.1/125 


170.7/130 


265.2/211 



^Parentheses indicate single-parameter 90% confidence regions. 
l-In the 0.3-10.0 keV band. 



Table 4. The number densities and energies of thermal and non-thermal electrons 



Filament No. 1 2 3 4 5 6 Total 



Thermal electrons 

E.M. [cm-3] 7.0x10^3 1.0 x lO^" 1.3x10^3 2.5x10^3 1.3x10^3 8.4x10^3 3.1x105" 

Density (nJ)^[cm-3] .... 0.43 0.51 0.26 0.41 0.30 0.60 0.45 

Total number (TVj) 1.6x10^5-* 2.1x105"* 5.0x10^3 6.1x10^3 4.3x10^3 1.4x10^4 6.8x105* 

Energy (Ej)b[ergs] 9.2 X lO*" 1.2 x 10-*5 2.8 X iC* 3.5 X lO** 2.4 x iC" 7.9 x ID** 3.9 x 10*5 

Non-thermal electrons'^ 

Density (n^'r)^[cm-3] .. . 4.7x10"* 5.5x10"* 7.9x10"* 7.1x10"* 6.5x10"* 7.6x10"* 6.2x10"* 

Total number (N^'^) .... 1.8 X 10^1 2.2 X 10^1 1.5 x 10^1 1.1 x 10^1 9.4 x lO^" 1.8 x lO^i 9.3 x lO^i 

Energy (S^T) [ergs] .... 2.0 X 10** 2.4 X 10** 1.7 x 10** 1.2 X 10** 1.1 x 10** 2.0 x 10** 1.0 x 10*5 

Injection efficiency (T)f I.lxl0"3 1.1 x 10"3 3.0 x 10"3 1.7 x 10"3 2.2 x 10"3 1.3 x 10"3 1.4 x 10"3 



*We assumed that the depth of the emitting volume is 1 pc and that the filling factor = 1. 
Tntegration from iJ^in = 0.24 keV to Bmax = 37 TeV (see text). 



